Quantum dimer model with extensive ground-state entropy on the kagome lattice 
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We introduce a quantum dimer model on the kagome lattice with kinetic terms allowing from 3 
to 6 dimers to resonate around hexagons. Unlike models studied previously, the different resonance 
loops appears with different signs (given by the parity of the number of dimers involved) . These signs 
naturally appear when performing the lowest-order overlap expansion (Rokhsar and Kivelson 1988) 
of the Heisenberg model. We demonstrate that the quantum dimer model has an extensive zero- 
temperature entropy and has very short-range dimer-dimer correlations. We discuss the possible 
relevance of this new kind of quantum dimer liquid to the physics of the spin-i Heisenberg model 
on the kagome lattice. 

PACS numbers: 75.10.Jm 75.50.Ec 



I. INTRODUCTION 



Quantum frustrated Heisenberg antiferromagnets are 
fascinating systems which can display a vast variety of ex- 
otic phases and phenomena. Systems with strong quan- 
tum fluctuations where no magnetic long-range order de- 
velops down to zero temperature ("spin liquids" loosely 
speaking) are of particular interest because they do not 
have direct classical analogs and are strongly interact- 
ing problems which resist to many simple theoretical ap- 
proaches. Focusing on two-dimensions and spin-i, two 
kinds of magnetically-disordered phases are well under- 
stood: valence-bond crystal (VBC) and short-range res- 
onating valence-bond (RVB) liquids. Both are character- 
ized by short-ranged spin-spin correlations but the VBC 
has long-ranged singlet-singlet correlations and gapped 
spin one excitations whereas RVB liquid has short-ranged 
singlet-singlet correlations, topological order and spin-i 
(spinon) excitations. 

Despite of intense theoretical ef- 

fortaWiWiWdVVVWV^ the physics of 

the spin-i kagomoi^ antifcrromagnctic Heisenberg 
model (KAFH) is still debated. For instance, there is 
still no consensus on the mechanisms which produces 
the unusually large density of singlet states that was 
observed numericallyiiiii 

Quantum dimer modelsiii£ (QDM) are effective ap- 
proaches to the phases of antiferromagnets which are 
dominated by short-range valence-bonds. These mod- 
els are defined in the Hilbcrt space of nearest-neighbor 
valence-bond (or dimer) coverings of the lattice and con- 
tain kinetic as well as potential energy terms for these 
dimers. These models can often be simpler than their 
spin parents and are amenable to several analytic treat- 
ments because of their close relations to classical dimer 
problems^, Ising models and 1*2 gauge theorjsi&SSiSI. 
These models can offer simple descriptions of VBC 17 
as well as RVB liquida22*2£ and a natural question is 
whether QDM can describe other phases, and in partic- 
ular whether they can describe phases with gapless sin- 
glet excitations. Motivated by the problem of the spin-^ 
KAFH, we investigated some QDM on the kagome lat- 



tice. Because of the corner-sharing triangle geometry, 
dimer coverings can be handled in a much simpler way 
(with pseudo-spin variable a 5 ! 20 ) than on other lattices. 
Exploiting this property we introduce a QDM (called //- 
model thereafter) with several interesting properties: i) 
The Hamiltonian allows dimers (from 3 to 6 at a time) 
to resonate around hexagons with amplitudes which have 
non-trivial signs. These signs are those arising when per- 
forming the lowest-order expansion (in the dimer over- 
lap parameter introduced by Rokhsar and Kivelson^ 7 -, 
see Sec. llVfl of the KAFH Hamiltonian in the valence- 
bond subspace^. These signs are the crucial difference 
with the solvable QDM we introduced previously^. For 
this reason also quantum Monte-Carlo simulations would 
face the well-known sign problem, ii) In addition to the 
topological degeneracy, a feature of dimer liquids, the 
ground-state has a degeneracy which is exponential with 
the number of sites, that is a extensive zero-temperature 
entropy, iii) The ground-states have short-ranged dimer- 
dimer correlations, they are dimer liquids. We studied 
the model through simple mean-field approximations as 
well as numerically and we propose a picture in which the 
system is critical (or at least close to a critical point). 

Because some parts of the paper are relatively inde- 
pendent, we will now summarize it so that readers may 
directly go to a specific part. In section [H] we review 
some results on the KAFH model. Although this pa- 
per is mostly devoted to a dimer model (sort of extreme 
quantum limit of the SU(2) spin-i model), we find it use- 
ful to review well established facts concerning the spin 
(Heisenberg) model and we motivate the QDM approach 
to the KAFH. In particular in Sec. Ill Fl we present nu- 
merical results (spectrum and specific heat) obtained by 
diagonalizing exactly (on finite-size systems) the Heisen- 
berg model restricted to nearest-neighbor valence-bonds 
subspacc. In Sec. IIIII we discuss general properties of 
dimer coverings on lattices made of corner-sharing tri- 
angles. These properties (existence of pseudospins vari- 
ables and their dual representation in terms of arrows) 
turn out to be useful to define and analyze QDM on 
these lattices, including kagome. In Sec. IIVI we explain 
the Rokhsar-Kivelson overlap expansion when applied to 
the family of lattices mentioned above. At lowest or- 
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der the kinetic energy terms have signs depending on the 
parity of the number of dimers involved. Ignoring the 
amplitudes and keeping only these sign (Sec. El, we get 
kinetic [ie non-diagonal) dimer operators n(h) defined on 
every hexagon h of kagome and which realize an original 
algebra: i) [i(h) 2 = 1, ii) they anticommute on neigh- 
boring hexagon and iii) commute otherwise. The rest of 
the paper is devoted to the analysis of the Hamiltonian 
defined as the sum of all the n(h). In Sec. IV Dl we start 
by solving exactly the dimer model on a one-dimensional 
lattice. It sustains critical (algebraic) correlations and 
has an extensive zero-temperature entropy. Although we 
did not succeed to find an exact solution to it, we were 
able to show (Sec. IV E|) that the kagome [i model also 
have such a zero temperature entropy. Some mean-field 
treatments are discussed in Sec. IV E 31 and a competing 
crystal-like phase is identified. In Sec. lVETl we introduce 
fcrmionic variables dual to the /i operators, in which the 
residual entropy is quite transparent. This formulation 
is reminiscent to that of the Z2 gauge theory in Ref. |2(]. 
The last section (IV1|) is devoted to numerical calculations 
on the kagome /i-model. 



II. SOME RESULTS ON THE KAGOME 
HEISENBERG ANTIFERROMAGNET 

In this section we review a few results concerning the 
Heisenberg model on the kagome lattice. 



A. Classical degeneracy 

The classical kagome antifcrromagnet attracted in- 
terest because of its unusual low-temperature proper- 
ties. These properties are related to the existence of 
a local and continuous degeneracy. Indeed, any spin 
configuration which as a vanishing total magnetization 
Si + 5*2 + S3 = on every triangle minimizes the Heisen- 
berg energy. Counting planar ground-states amounts to 
find the number of ways one can put A, B and C on the 
lattice so that each triangle has spins along the three dif- 
ferent orientations. This already represents an extensive 
cntropy^iS In a given planar ground-state one can look 
for closed loops of type A — B — A — B — Because on 
kagome such a loop has only C-type neighbors, rotating 
the spins of this loop around the C axis costs no energy 
and gives new (non-planar) ground-states. Chalker et 
alm& showed that all ground-states can be obtained by re- 
peated introduction of such distortions into the different 
parent planar states^ At low temperature this classical 
spin system has no magnetic long-range order (LRO) but 
exhibits diverging nematic correlations when the temper- 
ature goes to zero: although spin-spin correlations are 
short-ranged the planes defined by the three spins of a tri- 
angle are correlated at long distances. This phenomenon 
is a manifestation of "order by disorder" : thermal fluc- 
tuations selects ground-states with the largest number of 



soft modes and these are the planar ground-states. 



B. Absence of Neel long-range order 

It has been known for some time that the spin-i 
kagome Heisenberg antifcrromagnet has no Neel long- 
range order (LRO) at zero temperature. Early spin-wave 
calculations by Zeng and Elses^ indicated that magnetic 
order disappears when going from the triangular antifcr- 
romagnet to the kagome model. This was supported by 
numerical calculations of spin-spin correlations in finite 
kagome clusters up to 21 sites^. Two years later Singh 
and Husa^ performed a series expansion about an Ising 
limit and came to the same conclusion about the absence 
of magnetic LRO. 

Although the classical model has no Neel LRO at 
T = 0, the absence of such order in the spin-i case 
is not completely trivial because quantum fluctuations 
could select a particular type of ground-state. Sachde\*& 
showed in the context of a large N expansion that for a 
large enough value of the "spin" , a Neel LRO sets in (the 
so called y/3 x V3 structure). 

In 1993 Leung and Elser— pushed exact diagonaliza- 
tions to a 36 spins and confirmed the absence of Neel 
LRO. They also studied four-spin correlations (dimer- 
dimer) to investigate the issue of a possible valence-bond 
crystal (or spin-Pcicrls, or bond-ordered) phase made of 
resonating hexagons (see Fig.^). They found very weak 
correlations and suggested the existence of a liquid phase 
(they could not, however, definitely rule out the possibil- 
ity of a very weak crystalline LRO order). Nakamura 
and Miyashita— did Monte-Carlo simulations including 
N = 36 and N = 72 spins showed no kind of spin nor 
dimer ordering down to T ~ but found a low- 

temperature peak in the specific heat. 

On the analytical side, Sachde^ generalized the SU (2) 
model to an Sp(2iV) symmetry and worked out a large- 
N approach based on bosonic representations. He found 
quantum ordered phase with no broken symmetries and 
unconfined bosonic spinons. However this result does not 
directly explain^i the huge density of low-energy singlet 
states that was observed numerically and that we discuss 
below. 



C. Low-energy singlet states 

Lecheminant et ali& and Waldtmann et almA calcu- 
lated a large number of low-energy eigenstates for finite 
kagome clusters up to 36 sites. These results pointed to 
a large "residual" entropy at low temperatures. From 
their data the residual entropy per site can be estimated 
to be so ~ 0.21og(2). This number was obtained by 
counting the number of eigenstates in a finite (and non- 
extensive) energy window above the ground-state. This 
number was found to scale as ~ a N with a ~ 1.15. The 
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FIG. 1: Crystal of resonating hexagons (marked with 0) on 
the kagome lattice. Labels n = 0, 2, 4, 6 correspond to the 
possible resonance loop L n around each hexagon according to 
Table □ 

width of this energy window is expected to modify nu- 
merical prefactors^ but not a which is directly related 
to the an entropy per site sq = log (a). 



D. Specific heat 

The entropy change between zero and infinite tem- 
perature can be extracted from the specific heat c v (T). 
The first high-temperature (HT) series expansion for the 
kagome antiferromagnet was carried out by Elstner and 
YoungiS. This approach showed a huge entropy deficit of 
about 40% : / °° c v (T)/TdT ~ 0.61n(2). However, this 
direct evaluation of the specific heat from HT series is 
not accurate at low temperatures and they concluded the 
possibility of a low-temperature peak in the specific heat. 
A quantum Monte-Carlo simulation by Nakamura and 
Miyashitaai also found a low-temperature peak. Such a 
peak was also found by a decimation calculation^. A 
recent exact diagonalization work by Sindzingre el alSL 
also found such a peak in a 36-sites sample. An improved 
method of calculation of c„(T) from high-temperature 
series expansion, which is quantitatively accurate down 
to zero temperature in most frustrated magnets^ shows 
that about 20% of the total entropy is still missing at 
very low temperatures^, in agreement with exact diago- 
nalizations data. 



E. About residual entropies 

If a system has a number of states growing exponen- 
tially (with the system size) in a non-extensive energy 
window above the ground-state, it as a extensive resid- 
ual entropy at zero temperature. In such a case, although 
the ground-state can be unique on finite systems, it is 
in fact exponentially degenerate in the thermodynamic 
sensed. One can construct some simple models with an 
extensive residual entropy (the Ising antiferromagnet on 



the triangular lattice for instance) but it is usually lifted 
by almost any infinitesimally small perturbation. An ex- 
tensive entropy at T = is not a generic situation but 
instead requires some fine tuning (to zero) of all these 
perturbations. For these reasons we think it is unlikely 
that the spin-i kagome Hcisenbcrg antiferromagnet has 
a T = residual entropy. Consider some Hamiltonian 
7i(A) = H.q + XHi where TLq has an exponential ground- 
state degeneracy which is lifted by Hi- At small A the 
specific heat may have low temperature peak which en- 
tropy corresponds to the ground-state degeneracy of Hq. 
Upon taking the A — > limit, the temperature of the 
peak goes to zero as well. This is the picture we have 
in mind for the spin-i kagome antiferromagnet and this 
paper discuss a possible scenario in which the role of Hq 
is played by a quantum dimer model (defined as the /i- 
model in section W\ . 



F. Resonating valence-bond subspace 

Quantum dimer models can provide effective descrip- 
tions of some magnetically disordered phases of antifcr- 
romagnets. We first wish to motivate the restriction of 
the spin Hilbert space to the first-neighbor resonating 
valence-bond (RVB) subspace which has been used in a 
number of works^SiiSiffi for the kagome problem. This 
space is generated by all valence-bond states where spins 
are paired into first-neighbor singlets (dimers or valence- 
bond). Because spin-spin correlations are very short- 
ranged, it is rather natural to consider the ground-state 
wave function as a linear superposition of valence-bond 
states. The crucial point is to understand whether va- 
lence bonds beyond first neighbors should be included 
or not in the Hilbert space to (qualitatively) get the 
right physics. We will only partially address it in this 
paper. This first-neighbor RVB limit is the simplest sub- 
space which has an exponential number of states that 
could explain the proliferation of low-energy singlets ob- 
served numerically. In addition, this subspace provides a 
reasonably good variational energies. Zeng and Elscr 6 
and Mambrini and Mila™ computed the ground-state 
energy of the Hcisenbcrg Hamiltonian restricted to the 
first-neighbor RVB subspaceSi. For a sample of 36 sites 
their result (2 < Si-Sj >= -0.4218) is 3.8% higher 
than the exact ground-state energy obtained*^ in the 
full spin Hilbert space (2 < S l -S J > = -0.4384). Zeng 
and Elser were able to improve significantly this vari- 
ational estimate by a simple optimization of the dimer 
wave function in the vicinity of each defect triangle, but 
without changing the dimension of the Hilbert space. To 
our knowledge this "optimized dimerizations basis" is the 
best variational one for the kagome problem. It is also 
worth saying that in the fermionic large- TV extension of 
the Heisenberg model, first neighbor valence-bond states 
arise as degenerate ground-states in the N —> oo limit^i 
l/N corrections will then introduce a dynamics among 
these dimerized states. Marston and Zeng£ used such a 
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fermionic SU(N) extension of the Heisenberg model on 
the kagome lattice and found that such 1/N corrections 
could favor the crystal of resonating hexagons mentioned 
above (Fig. [JJ. 

A last argument for the first-neighbor RVB approach 
to the kagome problem is the fact that the spectra of 
the Heisenberg model projected into this subspace re- 
produce a continuum of singlet states as in the case of 
spectra computed in the full spin Hilbert space. This 
was first noticed by Mambrini and Mila-S on samples up 
to 36 spins and we extended their study to samples up 
to 48 sites. Figure El shows the exponential number of 
low-energy states in a finite energy window [Eq , Eq + 5] 
above the ground-state. We analyzed this exponential 
proliferation of energy levels as a function of the system 
size and as a function of the energy window. Although 
we have 7 complete spectra up to JV = 48 sites the de- 
pendence on the width of the energy window makes it 
difficult to give an precise estimate of the low tempera- 
ture entropy. For each value of 8 we plot the logarithm of 
the number of states in the window as in Fig. [21 A (least- 
square) fit is performed to extract the leading exponen- 
tial behavior when N — > oo. Error bars are obtained in a 
standard wayS In principle this procedure measures the 
zero-temperature entropy provided that 8/ J <C N. The 
result is summarized Fig. |3] Unfortunately, one cannot 
use too small values of 8 because discretization effects 
scatter the data when 8 is of the same order as the typical 
level spacing in the smallest samples. This is the reason 
for the increasingly large error bars we obtain when <5 is 
below 0.6 ~ 0.8 J (Fig. OJ. However, from these results it 
appears likely that a significant part of the total entropy 
is present at temperature much lower than the energy 
scale J. Indeed, the values of sq compatible with the set 
of data displayed Fig. [3 is 0.11og(2) < s < 0.21og(2). 
Only these values are within the error bars of all esti- 
mates from 8 = 0.4J to 8 = 1.2 J. 

Computing the specific heat is another method to 
look for a possible residual entropy. In the case of 
the kagome antifcrromagnct diagonalizcd in the full spin 
Hilbert space a low-temperature peak was observed^!, as 
well as in some experiment on a spin-3/2 kagome com- 
pound^. From its low sensitivity to an applied magnetic 
field, this peak was attributed (mostly) to non-magnetic 
singlet states. In this work we computed the specific 
heat of the kagome antiferromagnet in the first-neighbor 
RVB subspace. This calculation is done from the spec- 
tra obtained by numerical diagonalizations up to N = 48 
sites. The results shown are shown in Fig. 01 The max- 
imum of CV(T) around T = 0.7 J is almost converged 
to its thermodynamic limit. It corresponds to the onset 
of short-range correlations. For all sample sizes a large 
low-temperature peak is present at or below T = 0.07J. 
It is still size-dependent but its entropy roughly corre- 
sponds one half of the total entropy of the model (the 
total entropy per site of the RVB space is ^log(2)), in 
agreement with results of Fig. 01 The similarity between 
these results and those obtained in the full spin Hilbert 
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FIG. 2: Kagome Heisenberg antiferromagnet diagonalized in 
the first-neighbor RVB space. An exponential number of 
eigenstate is observed in the energy window [Eo, Eo+S] , where 
Eo is the ground-state energy. The results for two values of 
8 are shown. The full lines are quadratic least-square fits 
to the data for N > 18. The definition of J is such that 
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space is another support to the RVB approach. 

To summarize, we have reviewed several arguments in- 
dicating that the unusual low-temperature peak in the 
specific heat of the kagome antiferromagnet might be 
explained within the framework of a RVB space. We 
would like to conclude this section by mentioning that 
the spin-i Heisenberg model may have a large number of 
low-energy singlet states on other lattices made of corner- 
sharing triangles. This was observed numerically^ on the 
frustrated three-leg ladder shown Fig. 03 We wm come 
back to that model in Sec. IV Dl In Ref. l34l the squagome 
lattice was introduced and some low-energy states rem- 
iniscent of the kagome ones were identified in a large- V 
approach. A decimation method applied to this lattice 
and also predicts a low-temperature peak in the specific 
heat of the model^ 5 .. A numerical diagonalization study 
of the Heisenberg antiferromagnet on the Sierpinski gas- 
ket^& found a low-temperature peak in the specific heat 
as well. 



III. DIMER COVERINGS ON LATTICES MADE 
OF CORNER-SHARING TRIANGLES 

Before studying the restriction of the Heisenberg spin 
model to the valence-bond subspace we will introduce 
some properties of dimcr coverings on lattices made of 
corner-sharing triangles (including kagome) . A very use- 
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FIG. 3: For each width S of the energy window we fit the 
exponential increase of the number of energy levels to estimate 
the zero-temperature entropy so (as in Fig. Error bars 
come from the uncertainty of the least-square fits. 
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FIG. 4: Specific heat per site of the kagome Heisenberg model 
restricted to the first-neighbor RVB subspace. The lowest 
peak at T/J ~ 3.10~ 3 (TV = 48) is a finite-size artifact. 



ful property discovered by Elser and Zeng^ is that dimcr 
coverings on these lattices can be put in one-to-one corre- 
spondence with configurations of arrow variables. Also, 
this representation is intimately connected to the exis- 
tence of (Ising like) pseudo-spin variables^ 

The correspondence between dimer coverings on the 
kagome lattice and sets of arrows is illustrated in Fig. [3] 
Each arrow has two possible directions: it must point to- 
ward the interior of one of the two neighboring triangles. 
In a triangle, a dimcr connect two sites where the arrows 




FIG. 5: Arrow representation of a dimer coverings. 



point inwards. In a defect triangle (without any dimer, 
marked with * in Fig. |3J), the three arrows point out- 
wards. Therefore, at each triangle there is a constraint 
imposing that the number of incoming arrows is even. 

Dimer moves translate very simply in the arrow repre- 
sentation. One can easily verify that cr x (h) (see Ref. l2fl 
and appendix [SJ does nothing but flipping the 6 arrows 
sitting around hexagon h and that such an operation con- 
serves the constraint for all triangles. Any dimer moveSi 
is a product a x (h\)o x (I12) • ■ ■ where h\,h2,--- are the 
hexagons enclosed in the loops. This operation succes- 
sively flips all the arrows around h\,h2,- ■ •• The result 
does not depend on the order in which hexagons are 
flipped so the cr x operators obviously commute in this 
language. 



A. Medial lattice construction 

The arrow representation (as well as the pscudospin 
operators a x and a z introduced by Zcng and Elser^i - 
see appendix^) can be generalized to all lattices made by 
corner-sharing triangles. The kagome case is the simplest 
example in two dimension, another being the squagome 
latticed (Fig.|S}. The Sierpinski gasket — is an example 
of fractal structure of dimension between one and two 
also made by corner-sharing triangles. 

a) Start with a trivalcnt lattice H, that is a lattice 
where each site has three neighbors (full lines in 
Figs. EGO El and HU| The hexagonal lattice (Fig.0) 
is the simplest two-dimensional example. 

b) Construct its medial lattice K: Sites of K are, by 
definition, the centers of the links of H . The sites 
sitting on the three links of H connected to the 
same site of H are connected together. The medial 
lattice of the hexagonal lattice is the kagome lattice. 
Since H is trivalent, K is made of corner-sharing 
triangles. 
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FIG. 6: Medial lattice construction. Starting from a trivalent 
lattice (full lines) we construct a lattice which sites (black 
dots) are centers of the links. The sites of this new lattice are 
linked together (dashed lines) to form triangular plaquettes. 




FIG. 7: Kagome lattice (dashed lines and black dots) con- 
structed as the medial lattice of the hexagonal (full lines) 
lattice. The location of the pseudospins are indicated by up 
spins. 



c) Associate a pseudospin to each plaquettc of H (ie 
to hexagons of the kagome lattice in the example). 

In the following, we will use N for the number of sites 
in K, which is equal to the number of links in H . The 
number of sites in H will be 2N/3, which is equal to the 
number of triangles in K. The number of plaquette (or 
faces) in H is equal to the number of pseudospins, we 
write it -/V ps . For two-dimensional cases we can apply 
Euler's relation to the lattice H and we find (2N/3) — 
N + N ps =2 — 2g where g is its genus (g = 1 for a torus, 
g = for a sphere). 

ft ft J;' 

FIG. 8: Squagome lattice (dashed lines) as the medial lattice 
of the octagonal lattice. 




FIG. 9: A frustrated three-spin ladder (dashed lines)obtaincd 
as the medial lattice of (trivalent) two-leg ladders. 




FIG. 10: Another example of chain. 

B. Counting dimer coverings with arrows 

The number of dimer coverings of any lattice of type 
K (including, for instance, the one-dimensional examples 
of Figs. El and UUI) is 

My ■ - 2 N/3+1 (1) 

•I v dim. coverings V / 

This result can be obtained with the arrow representa- 
tion. Each arrow has two possible directions, which give 
2 N states. The fact that there can only be or 2 incom- 
ing arrows for each triangle introduces one constraint per 
triangle. There arc 2N/3 triangles but only 2^/3 — 1 con- 
straints are really independent, which gives Eq. ^ 

The fact that 2^/3—1 constraints are independent can 
be checked with the following argument. We focus on the 
trivalent lattice H on the bond of which the arrows live. 
First transform H into a tree by recursively cutting every 
bond that does not disconnect the lattice into two parts. 
The final tree is still trivalent so the number of leafs L is 
related to the number of vertices V by L = V + 2. Each 
bond gives two leafs when it is cut so that L/2 is the 
number of cuts. One can now set the arrows directions on 
the leafs. There are L/2 such independent arrows. Using 
the constraints associated to each vertex, the arrows are 
then determined on all the other bonds of the tree by 
progressively going from the leafs toward the root. It 
is simple to check that the last constraint encountered 
when reaching the root is automatically satisfied. From 
this we obtain 2 L / 2 = 2 v / 2+1 dimer configurations, which 
is equivalent to Eq. ^ since 3V = 2N. 

IV. FROM SPINS TO DIMERS: OVERLAP 
EXPANSION 

When restricted to the RVB subspace, the Heisen- 
berg model induces a complicated dynamics on valence- 
bonds. This dynamics is intimately related to the non- 
orthogonality of these valence-bond states that we de- 
scribe below. 
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A. Scalar product and loops 

The scalar product of two valcncc-bond states can be 
computed from their transition graph^! (loop covering 
obtained by drawing both dimerizations on the top of 
each other). We first need a sign convention for valence- 
bond states. A simple choice is to orient all the bonds so 
that all hexagons are clockwise^ (see Fig. lll|) . With this 
choice, the scalar product of two valence-bond states \a) 
and 1 6) is 




(a\b) = [] [(1/5 

Loops 



.L/2-1 



(-1) 



l+JVhex+L/2 



(2) 



where the product runs over non-trivial (of length > 2) 
loops in the transition graph of a and b, L is the length 
of the loop and Ah cx is the number of hexagons enclosed 
by a loop. For instance, the loop displayed Fig. ^2 h as 
A hex = 3. The factor (l/2) i/2_1 in Eq.His valid on any 



lattice whereas the signs (— 1) 



are associated 



to the corner-sharing triangle geometry^. The sign of 

L /2 

(a\b) would just be (—1) if all the bonds dimcrs were 
oriented clockwise around the loops. Consider the trian- 
gles on which a loop passes. We classify these triangles in 
four types (a, b, c and d) as follows. Some have two edges 
on the loop (types a and c in Fig. and the others have 
only one (types b and d). Some are inside the loop (types 
a and d) and the others are outside (b and c). Triangles 
of type a give two anticlockwise bonds and the (—1) fac- 
tors cancel. Every triangle of type b (resp. c) give one 
(resp. two) clockwise bond(s) and do not contribute to 
the sign. Each triangle of type d gives an anticlockwise 
bond and contributes by a factor (—1) to the scalar prod- 
uct. Let Nd be the number of such triangles. So far we 
have shown that sign [(a\b)] = H Loops {-l) Nd+L/2 . Using 
the arrow representation (and the associated constraint) 
one can show that (— 1) d = (—l) Nhe * +1 . The argument 
- not reproduced here for brevity - uses the fact the par- 
ity of arrows coming out and in of a given cluster of sites 
is fixed by the number of links and sites in that cluster, 
an thus to the number of hexagons enclosed by the loop. 

Rokhsar and Kivelsonii generalized the scalar product 
given by Eq.|3by giving two arbitrary fugacities x and y 
that couple to the number and to the length of the loops: 



fi 



x a.b 



n[ 

Loops 



L/2 



(-1) 



l+N hcx +L/2 



(3) 



FIG. 11: Transition graph between two dimer coverings. The 
bonds are oriented so that all hexagons are clockwise. The 
loop passes through four types of triangles: a,b,c and d (see 
text). 



B. Rokhsar and Kivelson scheme 

When restricted to the first-neighbor RVB space, the 
Heisenbcrg Hamiltonian induces a dynamics on valence- 
bonds. These valence-bond states are not orthogonal so 
we have a generalized eigenvalue problem. Orthogonal- 
ized valence-bond states la) can be obtained : 



i) = J2(^ 1/2 U\b) 



(4) 



and the matrix to diagonalize is 

Hf b {x) = (a\n\b) x (5) 

= ^(n-^w^'lwi&'U^/V (6) 



a'b' 



where H is the Heisenberg Hamiltonian H = ^<«> &i ' 

Sj. From now we will only deal with orthogonalized 
valence-bond states |o) but we will drop the tilde for 
clarity. H cS {x = 1/2) was diagonalizcd numerically to 
obtain the results of Fig. [21 and 0] (see also Refs. and 
l3(i|h Ti c f[(x > 0) is non-local and many dimers can hop si- 
multaneously to quite different configurations. However, 
since x = ^ < 1, the tunneling probability for such events 
decrease exponentially with the loop length. Up to order 
x n , H eS (x) is local and only contains terms with < n 
dimers. Following Rokhsar and Kivclson's work on the 
square lattice^, Zeng and Elses^ considered the small-x 
expansion of TC on the kagome lattice up to order x 6 . 
Up to a constant we have : 

n ° S (*) = - X! ha ^) \ L <») (L a \ + h.c. + 0(x 6 ) (7) 



x and y can be considered as formal expansion param- 
eters (x=l/2y = 2in the physical spin-i case). The 
choice y = 2 is usually adopted in the literatur e^' 1 ? but 
we will keep y explicit so that other cases may be con- 
sidered. When x —> the overlap matrix £l x becomes 
diagonal. 



where the sum on h runs over all hexagons, the sum on 
a runs over all the loops enclosing a that hexagon. The 
tunneling amplitudes h a (x) are given in Tab. |U These 
result reduces to Zeng and Elser's& when y = 1/x = 2. 
Notice that other terms of order x 6 exist and involve 6- 
dimer moves around two hexagons. 
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Loop a 



r2 




L/2 



h a (x) 



-3yx 3 = -3/4 



yx i (4-y) = l/4 



yx 4 {A-y) = 1/4 



2/^(4 - y) = 1/4 



2/a; 5 (2y-5) = -l/16 



ya; 5 (2y-5) = -l/16 



y:r 5 (2y-5) = -l/16 



^ 6 (6-3y) = 



+ 



+ 



+ 



+ 



+ 



+ 



+ 



TABLE I: The 8 different classes of loops which can surround 
an hexagon. Including all possible symmetries we find 32 
possible loops. L/2 is the number of dimers involved and h a 
is the tunneling amplitude (at lowest order) associated with 
each loop in a small- a; expansion of the Heisenberg model in 
the RVB space (see text). The value for the physical case 
y — l/x = 2 is given. 



Eq.0can be obtained from the scalar product formula. 
The latter is valid for any lattice K made of corner- 
sharing triangles (see section IIII A|) and Eq. [7| can be 
generalized to these lattices. Now hexagons are replaced 
by plaquettes of the lattice H. Consider a loop a encir- 
cling a single plaquctte. It has a length L and encloses 
N t triangles. The the amplitude h a (x) for that dimer 
move is 



h a (x) 



-Vi-xfl"- (L - yN t ) + 0{x L m J (8) 



where L m in (= 6 for kagome) is the size of the plaquettes 
of H . Unlike the squared or the triangular lattice case^ 2 ", 
no sign convention for the dimer coverings can turn the 
signs of the amplitudes h a all equal. 



V. QUANTUM DIMER HAMILTONIAN 

In Ref. |53 we introduced a solvable model: 

Ho = -J2\ L c)( L c l \+h. C . = -J2^(h) (9) 

h,a h 

where the pseudo-spin operators o~ x (h) are the kinetic 
energy terms defined in Table [I] (see also appendix lA"|l . 
This model is obtained by setting h a = 1 in Eq. [7| We 
showed^ that Eq. is completely solvable and is the 
prototype of RVB dimer liquid. It has a unique ground- 
state (up to a topological degeneracy) and its elementary 
excitations are pairs of gapped Ising vortices (visons). 
Here, we search for a model which, as Eq.03 is amenable 
to an analytical treatment and which could capture some 
essential features of the spin model. The first step in 
treating the frustration inherent to the Heisenberg model 
on the kagome lattice is to introduce non trivial signs in 
the dimer resonance loops. 



A. Definition of fi(h) 

We choose to keep only the sign of the leading terms h a 
of the dimer overlap expansion for the Heisenberg model. 
This leads to the definition of an operator /j,(h) at each 
pscudospin location h (hexagon centers in the kagome 
case): 



ft{h) 



32 

E 

a=l 



e(a) \L a ) (L a \ + h.c. 



where 



e(a) = (_1)Lcngth(a)/2 



(10) 



(11) 



and Length(a) is the length of the loop a. The ac- 
tion of /j(h) only differ from the pseudospin flip a x (h) 
by a sign: fi(h) \D) — ±a x (h) \D). This sign depends 
on the length of the admissible loop at h in state 
as indicated in Table [I] With this definition, if |c) 
is a (non-orthogonalized) dimer configuration, we have 
(c\ /i(/i) |c) > for any hexagon h. This can be used to 
define the signs matrix elements of fi independently of 
the orientation of the dimers. 

As for o~ x , (i can be simply expressed in terms of the ar- 
row representation introduced in Sec. 11111 Let 1, 2, • • • , 6 
be the sites of hexagon h and 1', 2', • ■ • , 6' the sites of 
the "star" surrounding this hexagon (see Fig. I12fl . The 
length of the admissible loop ao around h is related to 
the state of the arrows on 1', 2', ■ ■ ■ , 6' in the following 
way. If the arrow at site i' is pointing toward h, it shares 
a dimer with a site j belonging to the hexagon and ao 
will pass through i' and j. Let n out be the number of 
outgoing arrows. From this it is clear that the length of 
that loop will be Lcngth(ao) = 12 — n out . To summa- 
rize, the operator fi(h) defined in Eq. I1UI flips the arrows 
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5' 4' 



FIG. 12: The position of incoming arrows on sites l',2',--- 
defines the possible dimer move (full and dashed links) around 
the hexagon. The arrows at sites 1, 2, • • • , 6 are omitted for 
clarity; they take opposite directions in the dimerizations pic- 
tured by full and dashed dimers respectively. The arrows at 
1', 2', ■ ■ • are unchanged during this dimer move. 

of sites 1, 2, ■ • • , 6 and multiplies the configuration by a 
sign**: 

e(h) = (-1)»-*W/ 2 (12) 

The operators fx can be explicitly written with the help 
of a z and a x operators. Finding such an expression is not 
completely obvious since the a z are non local and depend 
on the reference configuration whereas the fi are local 
and independent of any reference state. The expression 
is derived in appendix FBI 

B. Commutation rules 

1. Operator fi 

The fi operators have the remarkable property of 1) 
anticommuting when they operate on nearest-neighbor 
hexagons but 2) commuting otherwise. This property 
is not obvious when looking at Ea. 1 101 and is most easily 
demonstrated with the help of the arrow representations. 
As mentioned previously the effect of fi(Ji) on an arrow 
configuration is to flip the arrows around the hexagon h 
and multiply it by a sign e(h) = (-l) n ™'(' 1 )/ 2 . Therefore, 
the action of fi(A) and fi(B) commute up to a sign. If the 
hexagons A and B are not adjacent, the signs e(A) and 
e(B) are unaffected by the action of fi{B) and fi{A) re- 
spectively, and fJ-(A) and fi{B) commute. For two neigh- 
boring hexagons A and B, the action of fJ,(A) affects the 
sign e(B) and conversely. There are two types of arrow 
configuration shared by the neighboring hexagons A and 
B, as shown in Fig. 113b and b. Both configurations 
have an odd number of outgoing arrows among the four 
external links which will be flipped by the successive ac- 
tion of fi(A) and fi(B). One of the signs e(A) and e(B) 
will therefore be changed by the action of the neighboring 
/i, but not the other. This means that upon acting with 
fi(A) and fi(B), the sign of the final configuration depend 
on the order we applied the two operators, and we find 
that on any configuration fi(A)fi(B) = —fi(B)fi(A). 




FIG. 13: Two possible states of a pair of triangles. 

2. Operator fi 

Another choice for the signs of the tunneling ampli- 
tudes turns out to be very useful. Consider the fl opera- 
tors defined by 

32 

/2(/i) = ^£(a)|L a )(L Q |+h.c. (13) 

where the signs e(a) are given in the last column of Tab. [I] 
Contrary to e(cv), e(a) counts the parity of the number of 
outgoing arrows only on one-half of the sites of the star: 
l',3' and 5'. Using the arrow representation and similar 
arguments as above, one can show that 1) fl operators 
anticommute when acting on nearest-neighbor hexagons 
2) commute otherwise. Most interestingly the fi and fl 
realize two copies of the same algebra that commute with 
each other: 

Vh,ti [n(h),fi{h')}=0 (14) 

The fi (and ft) operators have simple commutation 
relations with the pseudo-spin operators <j z introduced 
by ZE. By definition the fi operators can be written 
/i(/i) = e(h)a x {h) where e(h) is diagonal operator in 
the dimer basis. Because a z operators are also diago- 
nal in that basis they commute with any e and we simply 
have fi(h)a z (h) = -a z (h)fj,(h) and [/j,(h),a z (h')] = for 
h 7^ h' (see appendix for the commutation rules of a x 
and cr 2 operators). The same result holds for the fl. In 
Ref. |2(3 we identified a z (h) as the operator which cre- 
ates (or destroy) an Ising vortex (vison) on hexagon h 
and a non-zero expectation value (a z ) was interpreted as 
signature of confinement and vison condensation. 

C. Hamiltonian 

The main results of this paper concern the following 
quantum dimer model 

h 

where fi obey the commutation rules described in the pre- 
vious section. We studied this model both numerically 
and analytically. Numerically we diagonalizcd it on sys- 
tems up to 144 kagome sites (48 hexagons). The results 



10 




FIG. 14: The eight possible dimer moves around a square 
plaquette of the frustrated ladder shown in Fig. [5] The signs 
of the corresponding amplitudes in the /i-model are indicated. 

are presented in section IVTl The most striking feature 
of the spectrum is that all energy levels have a huge ex- 
tensive degeneracy of order 2 Np3 ^ 2 , where, in the kagomc 
case, N ps = N/3 is the number of hexagons. The de- 
generacy of the spectrum has its origin in the existence 
of the set of operators jl(h), which commute with the /x, 
and therefore with 7i M . The spectrum is organized in 2 Nc 
sectors labeled by the eigenvalues of the N c independent 
commuting operators which can be constructed from the 
jl, as explained in Sec. lVE"5l The eigenvalues are identi- 
cal in all the 2 Nc sectors. Another interesting feature is 
the existence of quantities which commute both with (i 
and jl, called in the following integrals of motion. They 
are constructed from products of fi (or alternatively jl) on 
straight lines drawn on the triangular lattice of pseudo- 
spins, as explained in Sec. IV E 71 The spectrum depends 
on the values of these integrals of motion. 

Using the operators jl, the dimer-dimer correlations in 
the model can be shown (Sec. IV E 6(1 to be short ranged. 



D. One-dimensional /i- model 

Before discussing the kagome case in more detail, it is 
interesting to look at the ^-model on the one-dimensional 
lattice displayed Fig.El As on any lattice made of corner- 
sharing triangles, the /z-model can be defined there. The 
different dimer moves (and their signs e(a)) are displayed 
in Fig. El Because it is one-dimensional, it is solvable 
and we will show that its spectrum exactly maps to the 
spectrum of the transverse-field Ising chain at the critical 
field (quantum critical point). 



1. Transverse-field Ising chain 

The mapping to the transverse-field Ising chain can be 
realized through a representation of the algebra of the \x 
by some pseudo-spin operators t x and t z defined by: 

T 2n = M2n T 2n+1 = fan 

t l = n ^ 2 '-! t 2«+i = n ( ie ) 

l<n l<n 

Using the (anti)commutation relations of fi and jl it is 
easy to check that if and r* anticommutc. The [i (resp. 



jl) operators only involve the Ising pscudospins on even 
(resp. odd) sites: 

M2n = T% n jl2n = T 2n+1 
M2n+1 = T 2n T 2n+2 fan+1 = T 2n+l T 2n+3 ( 17 ) 

And the [i Hamiltonian on the chain is now simply 

W M = - E T 2« - E T 2>2 Z „ + 2 (18) 
n n 

which we recognize to be that of a ferromagnetic Ising 
chain in transverse field at its critical point™. This model 
can be solved by a Jordan- Wigner transformation (re- 
duces to free fermions). Because the Ising chain is at 
its critical point, the dimer spectrum is gapless and sup- 
port linearly dispersive excitations at small momentum. 
In addition, (fJ-ifij) correlations decay algebraically with 
distance. The exponential degeneracy of 7i M is now trans- 
parent: only one half of the Ising spins (located on even 
sites) appear to be coupled by the Hamiltonian. However 
this entropy has a subtle origin: To write the the ^-model 
only in terms of the N ps /2 "coupled" degrees of freedom 
(in order to get rid of the entropy) one has to use op- 
erators (t|„) which are non-local for the original dimcrs 
(see Eg. I16|) . On the other hand, there are local oper- 
ators (the jl themselves) that do not change the energy 
and which create localized zero-energy excitations. 

We can make a comparison with another quasi one- 
dimensional model with extensive degeneracy: the spin-i 
Heisenberg model defined on a chain of coupled tctrahc- 
dra™. In that model some Ising-like degrees of freedom 
X = ±1 (spin chirality) do not appear in the Hamiltonian 
of the low-energy sector and the model has an extensive 
zero-temperature entropy. This situation seems analo- 
gous to the /x-model: the r| n+1 play the role of the \. 
However the important difference is that the T2 n +i are 
non-local in terms of the original dimers whereas the \ 
are local in term of the original spins. 

We will see in the next sections that the ^-modcl on the 
kagome lattice also has this important property that the 
"coupled" Ising degrees of freedom are non-local in terms 
of dimers. We are not aware of any other interacting 
quantum model exhibiting such kind of extensive degen- 
eracies. Their might be however an analogy with other 
systems with localized excitations such as in Aharonov- 
Bohm cages^£ or flat-band systems in general. There, 
the extensive degeneracy is due to destructive interfer- 
ences that prevent excitations from hoping on the lattice 
and thus gaining kinetic energy by derealization. This 
stresses again the role of the signs in the fi operators. 

2. Order parameter 

The mapping the to the transverse field-Ising chain 
suggests to introduce a different coupling for the (i on 
odd and even sites 

n(T) = -r]r^2„-E^n+i (is) 

n n 
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T 2n T 2n+2 



(20) 



so that r is the strength of the transverse field. It is 
known that the ground-state has (rf„) ^ ( T in T in+2) cx ~ 
cept at the critical point T = 1. In other words, the fi 
have different expectation values on odd and even sites 
for r ^ 1. This state (I ^ 1) has long-ranged (ni^j) 
correlations and is a crystal in the (fi) variables with 
a diverging structure factor S(n) = Y^, n (~ •"•)" (MoMn)- 
It is interesting to remark that from the Ising point of 
view the natural order parameter is the "magnetization" 
(r|„), which is non-zero for T < 1 phase but vanishes for 
r > 1. This order parameter is non-local in terms of the 
dimers (see Eq. I Kill . However this does not mean that 
for r < 1 the dimer system spontaneously breaks some 
hidden Ising symmetry. There is no such Ising symmetry 
in the dimer problem and the spurious Z2 redundancy 
was introduced in the Ising representation: reversing all 
the r z gives in fact the same physical dimer state, as it 
can be seen from Eq. El I n the dimer language this is 
a consequence of the fact that reversing all the arrows 
twice the is proportional to the identity. 



3. Heisenberg model on a frustrated ladder 

The spin-i Heisenberg model on the three-spin lad- 
der shown in Fig. |5| has been studied by Waldtmann et 
alm&. They considered a J — J' model where the hori- 
zontal coupling is J and the diagonal one (around square 
plaqucttes) is J'. From their numerical results (exact di- 
agonalizations for TV < 30 spins and DMRG for TV < 120 
spins) it appears that the system may be critical in the 
region 0.5 < J'/J < 1.25 (vanishing spin gap) and that 
a spin gap opens for J'/J > 1.25. Quite interestingly 
they showed that at J = J' the specific heat is quanti- 
tatively very similar to that of the kagome antiferromag- 
net and exhibits a low-temperature peak. The finite-size 
spectra also show a large density of singlet states above 
the ground-state (although probably not exponential in 
TV). These similarities suggest that the corner-sharing 
triangle geometry is an important ingredient to generate 
a large amount of low-energy singlet excitations and it 
would be interesting to investigate the possible relation 
between this three-spin ladder and the ^-model. 



4- Exact spectrum via fermion representation 

In the mapping to the transverse field Ising chain have 
neglected subtleties associated with boundary conditions 
as well as constraints on the physical space of dimer cov- 
erings. Indeed, the later has dimension 2 Ar / 3+1 whereas 
we used a representation of dimension 2 Nps = 2 N ^ 3 . Wc 
shall now present the full solution of the /z-model on the 
chain, using fcrmionic variables. 

Since we are interested mainly in the spectrum, we will 
realize the algebra of \ij operators in a space which has 



the right dimension of the dimer space, that is 2 Nps+1 . 
There exist an exact mapping between the fermionic 
states and the dimer states, but we do not insist on it 
here. We suppose periodic boundary conditions and take 
TVp S even. 

Let us introduce a pair of creation/annihilation 



fermionic operators Cj , c\ at each site j = 1, 



ps • 

They are equivalent to a pair of Majorana fermions, 
7j = (cj + Cj), jj = -i(ct - Cj ), 7 f = 7j 2 = 1. Wc 
construct the operators \ij as follows 



so that 



E 

3=1 



JV,: 



l'j 



3 = 1 



7373+1 



(21) 



(22) 



It is straightforward to check that two operators fXi, Hj 
anticommute if they are neighbors and commute other- 
wise; the same is valid for the operators fy, while \Xi and 
Jij always commute. 

We have not yet specified the periodicity conditions 
on the operators jj. Let us introduce two extra kinetic 
operators, fj,(u) and which move the dimers around 
the two edges (or TV ps -gones) of the chain 



/"(«)= E e (a) (l^a) + h-c.) 



u, d 



where a runs over the 2 Nps 1 possible loops of even 
length around the edge and e(a) = (_i)L<=ngth(a)/2^ 

They mutually commute and anticommute with all fij, 
j = l,...TV ps . Their product ^{u)n{d) commute with 
the Hamiltonian and it is an integral of motion, tak- 
ing the values ±1. Other two integrals of motion arc 
the products of fij on the even, respectively odd sites, 
Mo = ■ ■ ■ Vn ps -i, Me = M2M4 ■ ■ ■ Mjv P3 - A careful 
analysis shows that fx fi e fx(u)fJ,(d) = (— i) 3A W 2 on any 
dimer state. This is constraint on physical states shows 
that the dimension of the Hilbcrt space is 2 Nps+1 . Sup- 
pose now we are in the sector with fj,(u)/J,(d) = 1. Then, 
MoMe = — ( — l) Arpa ^ 2 MiM2 ■ ■ ■ MJV ps , and using the defini- 
tion from equation (|21|l we obtain that 7 i 7 at „+i = — lj 
which implies antiperiodic boundary conditions on the 
Majorana fermions. In the second sector, fi(u)fi(d) = — 1, 
the Majorana fermions have periodic boundary condi- 
tions. 

We have now all that we need to solve the model l|22|l . 
First, all the spectrum is degenerate 2 Nps ^ 2 times, since 
all the Np S Majorana fermions 7,- commute with the 
Hamiltonian. Then, the Hamiltonian [221 can be diago- 
nalizcd after a Fourier transformation. 



H/x = 2y^sinfc 7-^:7^ 



(23) 



where the sum is over momenta k = 27rn/TV ps , with n 
being an integer in the periodic sector a half-odd integer 
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in the antipcriodic sector, < n < N ps . The operators "fk 
with < k < it (resp. — it < k < 0) could be interpreted 
as annihilation (resp. creation) operators. Note that, in 
order to ensure the right commutation relations between 
creators and annihilators, we defined the Fourier modes 
as 7fc = l/y / 2N ps J2 n e lfc ™7n- The zero modes k = 0, it 
need separate treatment, but they do not appear in the 
Hamiltonian. After normal ordering, 

Hf, = 4 I sinfc l 7-fc7fc - 2 E I sinfc l • ( 24 ) 

0<fc<7T Q<k<ir 

The vacuum energy is easily calculated in the two sectors, 



E { - ] - 

£j Q — 



N pa /2-l 



E 



( + ) 



-2 E 

N p3 /2-l 

-2 E 



sm 



ir(2n + 1) 



ps 



sm(ir/Np a 



2irn 

w 



ps 



tan(7r/A^ ps ) 



(25) 



2. Three- sublattice representation 

One can use a different representation for the //'s on a 
smaller subspace, therefore removing part of the degen- 
eracy. Consider a decomposition of the triangular lat- 
tice into three sublattices A, B and C. All the /2(c) with 
c G C commute with each other (as well as with all the 
fx). They can be simultaneously diagonalized so that we 
can consider an eigenstate l^o) of these /2(c): 

M(c) l^o) = »?(c) v(c) = ±1 (27) 

In addition we may define some operators s(h) by: 

s{h) = n{h)ix{h) (28) 

They are diagonal in the dimer basis and commute with 
each other. We can project |^o) onto the subspace where 
s(c) = 77(c) for any site c G C: 



Since Eq ' < E^' , the ground state of the dimer prob- 

-2/sin(7r/iVp S ). In the 



n>i) = n 



lem has energy Eq = E [ 
thermodynamic limit E/N i 



1 + 7y(c)s(c) 



cited state is at Eq + ^; the rest of the spectrum can be 
constructed by making particle-hole excitations over the 
two fermionic vacua, according to Eq. 1241 The numeri- 
cal spectra, obtained in the dimer representation, are in 
complete agreement with the ones constructed from l2*4l 
In the thermodynamic limit the gap vanishes and the ex- 
citation spectrum of the 7^ is linear at small momentum. 



E. Kagome case 

1. Degenerate representation 

As we did in the one-dimensional case, it is natural 
to represent the /1 operators with Ising variables. The 
simplest representation uses one Ising variable t z = ±1 
at each hexagon: 



-2/tt. The first ex- which, bv Eg. 1271 is simply 



|V>i 



n 



1 + /j(c) 



IV>o) 



(29) 



(30) 



Now we consider the states generated by the action 
of the /i(aGA) and fj,(b€B) on \ipi). A basis can be la- 
beled by 2N ps /3 Ising variables 
follows: 



±1 and = ±1 as 



1 w+o 



I /i(&)^ l+ ^ 

.b£B 



l^l) 



(31) 

As we will now show, this subspace is stable under the 
action of any \i. This is obvious concerning the operators 
H located on sublattices A and B. On these sites we 
may define (non-diagonal) t x operators which reverse the 
value of t z at the corresponding site. With this definition 
we have 



n 



n 



(26) 



li(a) 
Kb) 



(32) 
(33) 



where the three unit vectors are at 120 degrees and 
relate a site to three of its neighbors on the triangular 
lattice. It is easy to check that this representation in- 
deed realizes the [i (and ji) algebra. One can in fact ex- 
press the (j, in terms of ZE pseudospin operators a z and 
find similar (although more complicated) relations (see 
appendix [BJ) . Notice in particular that t z anticommutes 
with [i and fl on the same hexagon but commutes with 
all the others, exactly as the u z do. This representation 
has (approximately) the correct size of the Hilbert space 
(^. 2^) but it does not show how many Ising variables 
decouple in this model, that is how large is the degener- 
acy (entropy) in this model. 



The t z t z ,t z „ comes from the anticommutation of /j,(b) 
with its three neighbors (a, a' and a") belonging to A 
when acting on a state like Eq. |^ Now we act with a 
/i(c) on \t%,t§). Upon moving /i(c) to the right through 
the /j(a) and n(b), the state picks a sign t z t z ,t z „t^t^,t^„ 
where a, a 1 ',■ • ■ ,b" are the six neighbors of c. Then we 
use the fact that /i(c) = \ipi) so that we finally get 
the following representation: 



H(a) 

Kb) 
H{c) 



(34) 
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This representation is independent of the 77(c). As a con- 
sequence, the spectrum of any Hamiltonian made out of 
\x operators (like TL^) will be at least 2 Nps ^ 3 fold degen- 
erate. As we show later that this degeneracy is in fact 
much larger ~ 2 NpB ^ 2 but this already exhibits an exten- 
sive residual entropy at zero temperature. 

3. Mean-field approximation 

The representation of Eq. 1341 suggests a simple vari- 
ational (or mean-field) approximation in which the sys- 
tem is in a tensor product of single-site wave-functions 

I*) = (g)f ps/3 ® |*b))- All the sites a £ A are in 

the same state as well as all the b e B: 

- cos(0/2)|T)+sm(0/2)||) (35) 
= cos(0/2)|T)+sin(0/2)||) (36) 

We have two variational parameters 9 and <j> an d the 
expectation value of the energy per hexagon is 

(Ha) = sin(0) 

(fi b ) = cos(0) 3 sin(^) (37) 
(fj, c ) = cos(0) 3 cos(<^>) 3 

E/N ps = - i «/0 + <ji b ) + (/i c )) (38) 

Minimizing _E we get 

(ji a ) ~ 0.2979 ; (/j,) ~ 0.3104 ; ( Mc ) ~ 0.7091 (39) 
£/iVp S = -0.4391 (40) 

Because (/i a ) 7^ 7^ (fj, c ) this state breaks the trans- 
lation symmetry. It has some crystalline order with re- 
spect to the (fi) variables. Notice, however, that such 
a state is not a dimer crystal (see Eg. I50|) . since it can 
be chosen to have zero dimcr-dimer correlations beyond 
a few lattice spacings. From the numerical diagonaliza- 
tions we estimate the exact ground-state energy to be 
E/N ps ~ -0.44±0.02 (see sectionlVTland Fig.EBJ, which 
agrees with the variational result within error bars. It 
is interesting to compare these energies with the energy 
that one get with a single-site wave- function derived from 
the translation invariant representation of Eq. 1261 The 
later approximation gives E/N ps = —0.3248, which is 
significantly higher. Two others mean-field states can 
be considered from Eq. [SB] On can use three different 
single-spin states \^a), \^b) an d \^c) on three sublat- 
tices. Minimizing the energy with respect to the three 
related angles we get E/N ps = —1/3. The corresponding 
variational state simply has (fi a ) = 1, (fJ-b) = (aO = 0. 
Enlarging the unit cell does not help to lower the energy. 
Indeed, using four sublattices leads to an even worse en- 
ergy (E/N ps = -1/4, (fi a ) = 1 and (fi h ^ d ) = 0). The 
fact that the degenerate representation gives rather bad 
energies can be explained from the fact in such states the 
Ising degrees of freedom (t x or r z ) are completely uncor- 
rected on different sites: (tjTj) = (Tj) (tj). The situation 



is different in a representation like Eq. |MJ In that case 
non-trivial nearest-neighbor correlations are present even 
in simple tensor-product states such as the one we con- 
sidered. 

Although the good variational energy given Ea.l40ldoes 
not prove that the system indeed spontaneously break 
the translation symmetry, it indicates that the /t-modcl 
on the kagome lattice is not very far from such a phased. 
The numerical results presented in section I VII indeed 
show that, at least at short distances, < [i(x)fj,(y) > 
correlations match the three-sublattice pattern. 



4- Three-sublattice [i-model 

As for the one-dimensional model, we can generalize 
the kagome /x-model by letting the couplings be different 
on sublattices A, B and C. 

H = -A4 Ma - X B ^2 W> - A c X! ^ c ( 41 ) 

a£A beB ceC 

and we focus on A^ + As + Ac = 1 and A^, As, Ac > 0. 
We determined the ground-state of the model within the 
mean- field approximation of Eqs. 1371381 The result is 
schematically shown Fig. ll5l (a qualitatively similar phase 
diagram is obtained with the degenerate representation). 
We obtain three phases. When Xa dominates the ground- 
state has = tt/2, (t*) = and (r 6 z ) ^ 0. When X B is 
the strongest we have 4> = tt/2, (t z ) 7^ and (r^) = 0. 
And close to the Ac point we have (t*) 7^ and {t§) 7^ 
0. Along the transition lines (dashed line in Fig. I15fl an 
expectation value (r z ) jumps from to a finite value so 
that the transitions are first-order in this approximation. 

The mean-field prediction for the topology of the phase 
diagram appears to be quite plausible. If there are indeed 
three phases, then, by symmetry, the point of interest 
Xa = Xb = Ac is the ending point of three transition 
lines. It is not clear whether the discontinuous character 
of these transitions is an artifact of the mean-field ap- 
proximation or a robust property. If the transitions are 
really first order, then the //-model spontaneously breaks 
the translation invariance and realizes a crystal in the \x 
variables. It may also be that Xa — Xb = Ac is a crit- 
ical point. Although we have no definite conclusion on 
this issue, some of the numerical results (susceptibility) 
presented in Sec. IVII suggests a critical point. 



5. Non-degenerate representation 

The constructions above use at most N ps /3 commut- 
ing fj, operators but this is not he maximum number N c 
of mutually commuting operators made out of the \x. In 
fact there are of the order of iV c ~ N ps /2 Ising degrees 
of freedom that decouple, as it can be seen from the fol- 
lowing argument. Divide the lattice in four sublattices 
A, B, C, D as shown in the Fig.UHl The "spins" jl(h) on 
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FIG. 15: Mean-field phase diagram of the /i-model on the 
kagome lattice with three different couplings Aa, A_b and Ac 
on three sublattices. Because the mean-field approximation 
breaks the symmetry which exchanges A, B anc C, the tran- 
sition lines (dashed) do not precisely match the symmetry 
axis. Such symmetry should hold in the real system and was 
restored here for clarity. Transitions are first order. 
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FIG. 16: Dividing the triangular lattice in four sublattices 
A, B, C and D. The open dots represent the "spins" entering 
to a bow-tie operator centered on the sublattice B. The black 
dots and the stars correspond to the positions of spins and 
centers of the bow-tie operators respectively, forming a set of 
mutually commuting operators. 



the sublattice A mutually commute. In addition, one can 
consider the "bow-tie" operators T = /2(l)/i(2)//(3)//(4), 
centered on sites of the B sublattice and involving neigh- 
boring "spins" of the C and D sublattices. These N ps /A 
bow-tie operators commute mutually and with the op- 
erators fi{h) from the sublattice A, which gives us an 
ensemble of N ps /2 commuting operators^!. These con- 
served quantities (/2 and T) can be used to eliminate de 
degrees of freedom on two sublattices in a similar way as 
we did for the three-sublattice representation. However, 
now, one cannot avoid obtaining a non-local representa- 
tion of the n algebra in terms of the spins living on the 
sublattices C and D. The procedure is briefly explained 
in Appendix. |0 



6. Arrow and dimer correlations 

Using the jl and T introduced above, we can show that 
arrows and dimers are completely uncorrelated in any 



state provided it is an eigenvector of the jl(a) and T(b) 
where a £ A and b G B are hexagons in the vicinity of 
the arrows. The proof is given below. 

Let us consider a triangle (1,2,3) on the kagome lat- 
tice. On these three sites we may define an arrow opera- 
tor ai whose value is 1 if the corresponding arrow points 
toward the interior of the triangle and otherwise. With 
this definition the dimer occupation number n\ 2 on bond 
(12) is n\2 = a\a 2 . Now we assume that we have three 



operators Oi (i = 1,2,3) which satisfy: 

Vi 0,0, = a. l 6 l (42) 

Vz ^ j 6 % a 3 = (1 - aj )6i (43) 

O t |0) = ± |0) (44) 



These relations just mean that 0\ flips the arrows on 
sites 2 and 3 but dot not touch the arrow 1, etc. |0) is a 
ground-state of the model which is also an eigenvectors 
for the Oi. Notice that if 0\ and 2 exist, 3 = 0\0 2 
is a valid choice. Under these conditions, we show that 
the rii are uncorrelated with any other bond which is 
unaffected by the three Oi. As we will explain, the Oi 
will be realized as local combinations of jl and T. From 
Eas. 1421 and 1431 it is simple to check that 

ni 2 d 3 = 3 (1 - n 12 - n 2 3 - n 31 ) (45) 

(plus cyclic permutations). Now let X be any operator 
that commute with the three Oi (later we will chose X = 
1 or X = nki where (kl) is a remote bond). Using Eas. 1441 
and 1451 the correlation (X riij) becomes: 

(Xn 12 ) - (0|Xn 12 (6 3 ) 2 |0) (46) 
= (0| 6 3 X (1 - ni2 - n 23 - n 31 ) 3 |0) (47) 
= (X) - (Xn 12 ) - {Xn 23 ) - {Xn 31 ) (48) 

Using the relations obtained by cyclic permutations and 
solving the three linear equations we find: 

(Xny) = \ (X) (49) 

Using X — 1 we get that the dimer density is j and using 
AT = riki (a remote link) we find 

1 2 

{nkiriij) = — = (n) (50) 

In order to complete the demonstration we still have to 
show that the operators Oi can be constructed for any 
triangle. Let the state |0) be an eigenstate of the /i(oG A) 
and T(fee B) described in the previous section. Depend- 
ing on its position relative to the four sublattices A, B, C 
and D, a triangle (123) (on kagome) will demand slightly 
different constructions for its Oi operators. For brevity 
we will only consider the case of a triangle located be- 
tween an hexagon a £ A and an hexagon b G B. It gen- 
eralizes straightforwardly to the other cases. Let i = 1 
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be the common site of hexagons a and b, i = 2 G a and 
i = 3 G b. It can be checked that the following choice 
satisfies Eqs. and 133 : 

Oi = f (ft) 

6 2 = (l(a)f(b) (51) 

6 3 = /i(a) 

(52) 

The result given Eq. 1501 shows that dimer-dimer cor- 
relations are extremely short-ranged^ in a ground-state 
like 10), it is a dimer liquid. The same result was found 
in Ref. |2(J in a gapped dimer model. Here, because the 
ground-state manifold has a huge dimension, it is likely 
that some perturbations are able to select (out of the 
ground-state manifold of Ti.^) states with some dimer or- 
der. Even in such a case we expect dimer-dimer correla- 
tions to be very weak in the vicinity of TL^. 



7. Fermionic representation and integrals of motion 

A version of the arrow representation can be given 
in terms of fermions. The advantage of such a formu- 
lation is that signs are naturally included. Unlike the 
one-dimensional case, such the fermionic representation 
cannot be used anymore to solve the model, since the 
Hamiltonian H15fl cannot be reduced to a quadratic form 
in fermions. However, it is still useful in gaining some 
insight about the model, for example it helps under- 
standing the extensive degeneracy of the spectrum. In 
particular it provides an argument which shows that the 
degeneracy is ~ 2 N >"' 2 not only on the kagome lattice but 
on any lattice made of corner- sharing triangles. Another 
advantage of this formulation is to understand the origin 
of the integrals of motions associated to products of (i 
along straight lines. 

We associate to each vertex of the kagome lattice a 
fermionic occupation number or 1. For a given dimer 
configuration, the corresponding fcrmion configuration is 
given by the following rules: a defect triangle has either 
occupation numbers 111 or 000. In the other triangles, 
there is one dimer connecting the sites with equal fermion 
numbers, 00 or 11. There is a constraint on the parity 
of number of fermions on each triangle, alternating on 
adjacent triangles, for example, triangles pointing to the 
right in fig. [S] have odd numbers of fermions, and that 
pointing to the left have even number of fermions. It is 
not difficult to see that the arrow variables and the occu- 
pation number variables are essentially the same and that 
the constraints on them arc of the same nature. In par- 
ticular, the counting of the degrees of freedom is similar 
for arrows and fermions, the constraint for each triangle 
eliminating one spurious degrees of freedom. 

On each kagome site there is a pair of cre- 
ation/annihilation fermionic operators Cj,Cj. As ex- 
plained above, the dimer space is equivalent to the 



fermion Fock space, with constraints on the occupation 
number on each triangle. 

As in the case of the chain, we transform the fermions 
into a pair of Majorana fermions, -fj = (cj + Cj), jj = 

—i(Cj — Cj). The algebra of operators /J.(h) can be realized 
now by using only operators jj 

(j,(h) =i , 
3 

where the arrow means that the product is oriented (for 
example, it starts at the leftmost site of the hexagon and 
runs clockwise). With this representation the /7-modcl 
contains 6 fermionic creation or annihilation operators. 
Since two adjacent hexagons have one kagome site in 
common, the associated fi operators anticommute. Dis- 
tant [i's commute, since they are constructed from even 
number of fermions, and each operators (i squares to 1. 
Similarly, one can construct 

(7(h) = — « JJ Ti » 

i 

with the product running in the opposite direction to 
that of (i. These operators obey the same algebra as 
(i(h) and commute with the whole set of operators (i. 
The symmetric role played by the operators 7^ and 7, 
suggests that one half of the degrees of freedom are not 
affected by the action of the Hamiltonian Eq. ^| and 
therefore that the degeneracy of the spectrum is of the 
order of 2 JV *»A 

Let us note that the operators (i(h) and p,(h) leave the 
constraints on the occupation numbers invariant, since 
they change by two the occupation number on each tri- 
angle. More generally, products of 7 on loops which visit 
each triangle an even number of times also leave the con- 
straints invariant and their action could be translated 
in the dimer language. The most general operator which 
preserves the constraint (dimer space) can be constructed 
from products of 7 and 7 on loops visiting each triangle 
an even number of times. These operators are the equiv- 
alent of Wilson loops in a gauge theory. Triangles where 
the constraint is not obeyed can be constructed by action 
of strings of 7 operators (they could be useful to describe 
dimer configurations with defects, and therefore to intro- 
duce holons or spinons). And lastly, the vison creation 
operator is naturally constructed in this language as the 
product of 7^7^ on a string. 

Integrals of motion. Let us analyze the case of closed 
systems with the topology of the torus, which is the ge- 
ometry we used for the numerical diagonalization. The 
lattice on which the operators (i live is triangular and it 
is made by the centers of hexagons of the kagome lattice. 
Call toi, 771.2, "13 the number of closed lines in each el- 
ementary direction on the lattice, having length m, 772, 
n.3 respectively^, so that N ps = mini = 7712^2 = 771377,3. 
We denote by M the total number of such lines, M — 
777i + n r >>2 + ITI3. We label by L^, k = 1, . . . , M the closed 



16 



straight lines on the triangular lattice, and lk and lk+i 
the lines on the corresponding kagome lattice bordering 
the line Lk (with obvious periodicity conditions). Asso- 
ciated to each kagome line lk , we can define the following 
operators 

= n ^ » r fe = n tj > 

where the products run over the sites j of the kagome 
line Ik , with some ordering indicated by the arrow. Since 
every line lk visits any hexagon or 2 times, the operators 
Tk and I\. commute with all the operators fi{h) and fj-(h). 
From the correspondence with the dimer states, we know 
that Tfe and Tfc correspond to dimer moves on non-trivial 
loops around the torus, so they change the topological 
sector^ 

In each topological sector we can construct the follow- 
ing integrals of motion, associated to the closed straight 
lines Lk on the triangular lattice 

h = (-*)'Trrr +1 = n > 

heL k 

h = (-i)"fr +1 fr = n aw 

heL k 

where n denotes the number of hexagons on Lk and the 
arrows denote the ordering in Tk, IV The the third mem- 
ber of both equalities is invariant by circular permutation 
of the sites on the line Lk- 

The two sets of integrals of motion {Ik}, {Ik} are not 
independent. To check this, we use the commuting vari- 
ables 

s(h) = n(h)fj(h) , 

measuring the parity number of fermions around the 
hexagon h, 

s(/i)=n^7i = H) 6 II( 1 - 2n i)- 

3 3 

Then, 

ijk = (-i) n n s w = ■ 

where the last equality is a consequence of the constraint 
of fermions number around the two types of kagome tri- 
angles. Moreover, integrals of motion corresponding to 
the lines Lk with the same orientation on the lattice are 
not all independent, 

mi 

n h = Hr 1 " 1 = (-i) N " , (53) 
fc=i 

an similarly for the other two orientations. Such a con- 
straint on physical states is not unexpected, since the 



product of the lines with some orientation contains all 
the operators n(h) exactly once, that is all the operators 
7i twice, so it has to be proportional to the identity. 

In conclusion, there are at most M — 3 independent 
integrals of motion, where M is the number of different 
closed straight lines one can draw on the triangular lattice 
wrapped on the torus (in some cases, the M — 3 lines are 
not independent and some of the quantities Ik can be 
written as products of the others). 

Basis for the Hilbert space. The commuting quanti- 
ties s(h) and Ik can be used to label the states in the 
Hilbert space. Let us first check that we obtain the right 
dimension. Due to the M — 2 independent constraints 

f[hGi s (^) = 1> on ly ^p s ~ ~ 2) 01 t ne operators 
s(h) arc independent. Both type of operators s(h) and 
Ik can take only two values, ±1, or ±i for Ik on a line 
of odd length. Taking into account the topological de- 
generacy, the number of states in the Hilbert state is 

4 x 2 N P s-(M-2) x 2 M-3 = 2 Ar ps + l ^ which ig thc right 

dimension of the dimer space (in the case when there 
are extra relations between the lines Lk, the number of 
independent integrals of motion and the number of con- 
straints on s(h) arc diminished by the same number). 



VI. EXACT DIAGONALIZATIONS OF THE 
^-MODEL 

Because the non-diagonal matrix elements of the \i 
operator have different signs, Ti.^ is not appropriate for 
large-scale Monte-Carlo simulations and, instead, we per- 
formed some numerical diagonalizations of thc Hamilto- 
nian. For systems up to iVhex = N ps = 16 hexagons 
(N = 48 sites) we diagonalize it directly in the basis 
of dimer coverings (which dimension is 2 Nps+1 ) by using 
all lattice symmetries. For larger systems (N ps = 20, 
24, 28, 36) we diagonalize the Hamiltonian in a repre- 
sentation where the extensive degeneracy (due to the 
N c Ising quantities that commute with every fii) have 
been removed (see appendix Q . The diagonalization is 
performed separately for each sector defined by the con- 
served quantities I. For the largest system (N ps = 48, 
JV = 144) we use a Lanczos algorithm in this non- 
degenerate representation to obtain the first energies and 
wave-functions. Thanks to these numerous symmetries 
the largest vector size is only ~ 10 6 . The lattices we 
used are displayed Fig. IT7I 



A. Spectrum 

The ground-state energy per hexagon is plotted in 
Fig. El From this data we can estimate that < fj, >~ 
—0.44 ± 0.02 in the thermodynamic limit. It is interest- 
ing to compare this value with the energy of a simple 4-/i 
problem Ti = fix + \ii + fi3 + fi4 (with periodic boundary 
condition so that every site is neighbor of thc 3 others), 
which has < fi >= — i in its ground-state^ The energy 
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FIG. 17: Finite-size lattices used for the numerical diagonal- 
izations of the /x- model. Each dot is an hexagon of the kagome 



lattice. N v 



|iV is indicated. 



The lattices on the top are 



compatible with a three sublattice structure whereas the oth- 
ers (bottom) are not. All lattices except N pa = 10, 14, 16 , 16 c 
are compatible with the four-sublattice structure. 
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FIG. 18: Top: Ground-state energy per site. The dotted line 
is a least-square fit of the data for N ps > 16. Bottom: energy 
gap between the ground-state multiplet and the first excited 
state. The dotted line is a guide to the eye (~ e~ aNps ). When 
different values are plotted for the same N ps , they correspond 
to different shapes (see N ps = 16 and 24 in Fie. 1171 . 
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gap between the ground-state multiplet and the first ex- 
cited state is also shown in Fig. ^] This quantity prob- 
ably goes to zero in the thermodynamic limit. Fig. 1191 
shows the gap as a function of the ground-state energy 
per hexagon. It appears that the samples with the largest 
gap are those which energy is significantly lower that the 
thermodynamic estimate (< /U >~ —0.44 ±0.02). This 
also points to a gapless spectrum in the limit of large 
systems. 

The dispersion relation of the first excited states usu- 
ally provides some useful insight. However, due to the 
extensive degeneracy, this brings no information for /i 
model because the dispersion relation can be shown to 
be perfectly flat. Let |fc) be an eigenstate with mo- 
mentum k and energy Ek ■ The new state | k + q) = 
Ylh eH , rh fi(h) \k) has, by construction, momentum k + q. 
Because the fi commute with H^, \k + q) is still an eigen- 
state with energy Ek ■ This property is just a consequence 
of the fact that acting with a /2 creates a localized zero- 
energy excitations. 

From these data we propose two possible scenarios: 
1) An additional ground-state degeneracy associated to 
some spontaneous symmetry breaking in the thermody- 
namic limit. As we will explain, it may be that the 



FIG. 19: Gap as a function the energy per site for N ps > 12 
(data of Fig. 11811 . The line is a guide to the eye. 



ground-state orders in the three-sublattice pattern dis- 
cussed in section IV E 31 As we will explain, this scenario 
does not seem to be the most likely. In particular, if the 
system was a three-sublattice crystal in the fi variables, 
the spectrum would have a ground-state with a small 66 
quasi-degeneracy (ignoring the exponential factor com- 
ing from the jl degrees of freedom) separated by a gap 
to higher excitations. This does not seem to be the case 
since we could not identify a small set of energy levels 
adjacent to the ground-state above which a significant 
gap could exist. 2) In the second scenario the low-energy 
states may correspond to a gapless mode of excitations 
in the system. Although we have no precise theoretical 
prediction for the nature of such (critical) excitations in 
this 2D dimer model, gapless excitations are reminiscent 
of the one-dimensional analog (which has fcrmionic crit- 
ical excitations at low energy see section IVDl concerning 
the fj,- model on the chain). According to the single- mode 
approximation discussed in paragraph IVI Dl these gap- 
less excitations would only exist at a finite value of the 
momentum. 
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B. Correlations 



We looked for the possibility of long-ranged Hifij corre- 
lations in the ground-state. We define a static structure 
factor S(k) in the usual way: 



s(k) = ^<(W- fc |o) 



iy ps 



-ik-{ri —rj ) 



\(H fij |0) 



(54) 
(55) 




where 



(56) 



These calculations were done numerically in a reduced 
Hilbert space where iV c ~ N ps /2 conserved quantities 
(made out of ft and T operators) are fixed to be ±1 (non- 
degenerate representation). By construction the spec- 
trum does not depend on these choices (that is the ori- 
gin of the entropy) but it is also possible to check that 
(fii Hj) correlations do not depend either on the sector. 
However we stress that it is in principle possible to have 
different correlations in a ground-state which would be a 
linear combinations of the ground-states of different sec- 
tors. This is similar to the question of dimer-dimer cor- 
relations discussed previously. We have not investigated 
these effects which are related to the possible ordering 
pattern which may be selected by small perturbations in 
the ground-state manifold. 

The results are summarized in Figs. |21 IM1 and HH 
Fig. |2] clearly indicate that the most important cor- 
relations appear at the border of the Brillouin zone. 
More precisely the corners of the Brillouin zone ks = 
(±47r/3,0) and the middle points of the borders of the 
Brillouin zone fc^i = (0,27r/\/3), kA 2 = (it, — 7r/v3) 
and kA 3 = (tt, 7t/v3) are the reciprocal lattice points 
where the correlations are the strongest, fcs correspond 
to a 3-sublattice structure whereas fc J 4 i is related to a 

2- sublattice (stripe-like) order. A (weak) tendency to a 

3- sublatticc ordering can been seen directly in Fig. I2UI 
which represents real-space correlations in the ground- 
state of the N ps = 48 sample (144 kagome sites). Almost 
all the sites with a positive correlation (black circles) are 
locate on the same sublattice (according to a 3-sublattice 
decomposition) as the reference site. 

To check whether these correlations could remain long- 
ranged in the thermodynamic limit we plot the S(k)/N ps 
as a function of N ps (see Fig. I22fl. As a result, S(k)/N ps 
seems to extrapolate to a very small (possibly 0) in the 
limit of large systems. This suggest neither 2- nor 3- 
sublattice "crystalline" order in the expectation values 
of the fi operators. However, the data at k = ks should 
be compared with the mean-field state described in sec- 
tion According to the expectation values given by 
Eq.Elwe should have S{k)/N ps ~ 0.0182 in the thermo- 
dynamic limit. While the extrapolation of the numerical 



FIG. 20: Correlations {[ioHi} c = (/-toMi) — (^o) 2 in a 48- 
hexagon system. The site is at the center. The radius of 
the circles is proportional to | {^of^i) c \- Empty circles indicate 
negative correlations and the black ones are for {fj,ofJ.i) c > 
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FIG. 21: Structure factor S(k) = {fj,~kfJ>k) represented in 
the first Brillouin zone of the triangular lattice for N ps = 48 
(144 kagome sites). The radius of the circle is proportional to 
S(k). S(k) has a trivial divergence at k = which is due to 
the fact that the operators < fii > is non-zero at every site. 
This peak at k — is not represented here. 



results of Fig. [221 can n °t distinguish such a small order 
parameter from a disordered (or critical) phase, the pre- 
diction of the mean-field approximation (dashed line in 
Fig. I22|) turns out to be significantly different from the 
exact ones. In the mean- field approximation S(kB) is 
given by S(k B )/N ps ~ 0.0182 + 0.6279/A^ (the l/N ps 
contribution is just the local contribution of a given site 
and its six neighbors). On the other hand the exact 
value of S{kB)/N ps decays much faster with N ps . This 
means that the reduction of the crystal order parame- 
ter S(kB)/N ps with the system size is mainly caused by 
long wave-length fluctuations rather than by local con- 
tributions. This is a serious indication that the three- 
sublattice crystal is unstable with respect to these fluc- 
tuations. 
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FIG. 22: Top: structure factor at k — kA- Bottom: Struc- 
ture factor at k — ks- The dotted lines are obtained by a 
least-square fit of the form a + bNpJ~ + cN~ s 2 . The dashed 
line is the mean-field fEg. l4UH prediction for S(kB)/N pa . The 
quick reduction of the crystal order parameter with the sys- 
tem size (compared with the mean-field result) suggests that 
the crystal is unstable to fluctuations. The analysis of the re- 
lated susceptibilities confirms that S(k)/N ps is indeed likely 
to extrapolate to zero when N p3 — > oo. 



C. Static susceptibility 

To get more insight on the possibility of some crys- 
talline order in the \i variables we calculated the static 
susceptibilities x(fc): 



x(k) 



1 d < Hk + M-fe > 



2K 



ps 



dX 



(57) 



where A is the strength of a infinitesimal symmetry- 
breaking perturbation: 



7-i\ = ^2 to ~ \^ (to; + to-h) 



(58) 



x{k) is obtain numerically by measuring the expectation 
value of fXk + M-fc m the ground-state of the Hamilto- 
nian Eg. 1581 in the presence of a small perturbation. The 
susceptibility is obtained by extrapolating the result to 
A = 0. 

The static susceptibility is a rather sensitive probe 
since it must diverge as N 2 S in systems that sponta- 
neously break the translation symmetry in the thermo- 
dynamic limit 4i On the other hand it remains finite if 
there is no ordering at the corresponding wave-vector. 
The results for k = kA and k = fcs are displayed in 
Fig. x(k) shows no tendency to d iverge at k — kj\ 
and the increase with N ps of x(^b) is significantly slower 
than ~ N 2 , as suggested by the rather good fit obtained 
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FIG. 23: Static susceptibility x(k) for the wave vectors 
(k = kA and k = ks) where correlations are the strongest. 
Long-range order with spontaneous symmetry breaking would 
imply x ~ N 2 , which does not appear to be satisfied. The 
dotted line is a least square fit of the form x(^s) — aN + b. 



with x(^b) — aNp S + b (dotted line in Fig. 125)) . For these 
reasons we think that the system does not develop long- 
ranged correlations in the thermodynamic limit. 
The data for x(fce) (which neither diverges like N 2 S not 
stay constant) might be interpreted as a proximity to 
a critical point where a thrcc-sublatticc structure would 
appear. 



D. Long-wavelength fluctuations 

We now turn to the analysis of the long-wavelength 
fluctuations in the system. The structure factor S(k) is 
represented as a function of |fc| in Fig. [21 S(k) seems 
to vanish as least as S(q) ~ \k\ 2 and probably faster. 
A ~ |/c| 4 behavior looks plausible and is reminiscent of 
quantum Hall effect 42 As explained before, the disper- 
sion relation is flat in this model. However, one may be 
interested in the excitations that can be created by the 
action of the fx operators only (excluding the ft). Such a 
variational excited state with momentum k can be con- 
structed in the from the ground-state in spirit of the sin- 
gle mode approximation: 

\k) = ti k |0) (59) 

The energy (relative to the ground-state) of \k) is 

2WZ iM-k, [W,/ifc]] 



L0(k) 



S(k) 



(60) 



Since the numerator (oscillator strength) behaves, as 
usual, as ~ \k 2 \ at small \k\, we find that \k) is not a 
low energy excitations as soon as S(k) vanishes like \k\ 2 
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FIG. 24: Long-wavelength behavior of the structure factor 
S(k) for different sample sizes. 

or faster - which seems to be the case. This is suggestive 
of a non-zero gap for zero-momentum excitations. 



VII. DISCUSSION AND CONCLUSIONS 

We have introduced a QDM on the kagome lattice with 
a kinetic energy which allows from 3 to 6 dimers to res- 
onate around hexagons. The crucial difference with pre- 
vious QDM is that dimer move with amplitudes which 
have non-trivial signs inherited from the underlying spin- 
\ model. Exploiting the algebraic properties of the con- 
served quantities (Jl operators) we showed that the model 
as an extensive entropy at zero temperature - g ln(2) per 
kagome site - and is a dimer liquid. 

The starting point of this study was the spin- ^ Heisen- 
berg model on the kagome lattice. Concerning this prob- 
lem our main result is that a high density of singlet states 
at low energy might have a real quantum origin and may 
not just be the remainder of the local degeneracies of the 
classical model. The mechanism which produces the en- 
tropy of the fi model is new in the sense that one has 
to use non-local degrees of freedom in order to compute 
the spectrum in a representation which eliminates the 
degeneracy. We are not aware of any other model with 
a similar behavior. Going from the KAFH model to the 
fi model we made some crude approximations. The first 
one was to reduce the spin Hilbert space to a short-range 
RVB one. We provided several arguments to support 
this approximation but some additional studies would be 
required to analyze this question further. The second 
drastic approximation was to reduce the dimer dynam- 
ics induced by the Hcisenbcrg interaction to that of the 
/i-model and its signs. This can be qualitatively justified 
for the KAFH in the temperature regime corresponding 



to the low-temperature peak of the specific heat. If that 
picture is correct the degrees of freedom involved in that 
peak would correspond to the ~ N ps /2 = N/6 uncoupled 
degrees of freedom of the /i-model. The //-model can be 
defined on any lattice made of corner-sharing triangles, it 
could therefore provide a rather general explanation for 
a large entropy at low temperature in the correspond- 
ing frustrated spin-i models. Determining if some order 
eventually develops at much lower temperatures would 
amount to analyze a degenerate perturbation theory in 
the ground-state manifold of the //-model. 

An important question is to know whether the //-model 
realizes a new phase or if it is at a critical point. We have 
shown that the most serious candidate for an ordered 
phase, if any, is the three-sublattice crystal. However we 
gave several indications (spectrum, correlations and sus- 
ceptibility) suggesting that it is not stable. Instead we 
suggest that the system might be at a critical point. If we 
think of a dimer model as a system of hard-core bosons, it 
is interesting to compare our findings with some known 
bosonic phases. Let us first come back to the gapped 
RVB state realized in the solvable QDM of Ref.Hj. That 
state, which is the equal-amplitude superposition of all 
dimer configurations (Rokhsar-Kivelsoni£ state), is very 
similar to a Bose condensate in the sense that its wave- 
function can be obtained by putting all dimers in the 
same zero-momentum stated. The important difference 
with a conventional supcrfluid is that the dimer model 
has no t/(l) gauge symmetry (and therefore no gapless 
"sound" mode or conserved integer charge) but a discrete 
Z2 gauge symmetry^. With that comparison in mind the 
ground-state of the /i-model would neither be a conden- 
sate nor a crystal but it has gapless excitations. This is 
rather unusual in a model which has no continuous sym- 
metry at all. In addition our model has a structure factor 
S(q) for (/i/i) correlations which decays \q\ 2 or faster in 
the limit q — > 0. In a single-mode approximation this 
would imply a gap for q = excitations. 

In order to study the spectrum of the /i-model we used 
a representation in which the degrees of freedom respon- 
sible for the extensive entropy are frozen. The /i-modcl 
is local in the dimer variables but the effective Hamilto- 
nian describing the non-degenerate spectrum turned out 
to be non-local in terms of the original dimers. These ef- 
fective long-ranged interactions between physical degrees 
of freedom might be an important ingredient for the ap- 
pearance of an exotic liquid. From this point of view 
there might be some similarities between the /i-model and 
some two-dimensional quantum systems of Bosons with 
long-ranged interactions. Exotic liquid states which are 
not superfluid have been proposed for these systems, in- 
cluding a quantum hexatic phas o 44 i 45 . On the other hand 
a critical Bose fluid can exist without long-range interac- 
tions and such a phase can be stabilized by cyclic ring- 
exchanges A striking feature of these new phases 
is the existence of gapless excitations along lines in the 
Brillouin zone, as in a Fermi liquid. We have no direct 
indication of such a behavior in the kagome /i-model, ex- 
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ccpt for the existence of a fermionic representation, but 
such a scenario is certainly an interesting possibility that 
should be tested in future studies. 
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APPENDIX A: ZENG AND ELSER'S 
PSEUDOSPINS REPRESENTATION OF DIMER 
COVERINGS 

Zcng and Elser (ZE) realized*^ that a close correspon- 
dence between Ising configurations of pseudospins sitting 
on hexagons and dimer configurations on the kagome lat- 
tice (within a given topological sector) could be used. We 
used this representation in a previous workSS to define an 
exactly solvable quantum dimer model. In this section we 
review the pseudospin representation. 



1. a z component 

We first need an (arbitrary) reference dimer configu- 
ration | Z?o)- We will associate a pair of pseudospin con- 
figurations {<7^ = ±l} h(Ehcx to any dimer configuration 
\D) (belonging to the topological sector of \D )) in the 
following way: 

• Draw the loops of the transition graph of (Dq\D). 

• These loops must be considered as domain walls 
separating hexagons where the pseudospins are up 
and hexagons where they are down. This can be 
done in a consistent way because D and Dq arc 
supposed to be in the same sector and any closed 
path will necessarily cross an even number of do- 
main walls. 

• There is a two-fold redundancy in the above pre- 
scription because up and down hexagons can be ex- 
changed without changing the loop pattern. Since 
there is no natural way to decide where is the 
interior and where is the exterior of a closed 
loop on a finite sample, a pseudospin configu- 
ration {(j^ = &(h)} and its reversed counterpart 
{cr| = —a(h)} represent the same dimer covering. 

This establishes a one-to-one correspondence between 
dimer coverings of a given topological sector and pairs of 
pseudospin sates related by a global pseudospin flip. The 
proof can be done in two steps: 

• If two dimer configurations D and D' are associ- 
ated to the same pseudospin state (up to a global 
pseudospin reversal), they are identical. The tran- 
sition graph (D\D') can be viewed as the "differ- 
ence" between graphs (Dq\D) and (Dq\D'). This 



means that the (D\D') will have loops separating 
regions where the pseudospins coincide in D and D' 
and region where they are different. Since D and 
D 1 have the same a z on every hexagon, (D\D') can- 
not contain any hexagon and therefore contains no 
loop at all. 

• Any pseudospin configuration has a corresponding 
dimer state. The transition graph between Dq and 
the dimer configuration we are looking for will sep- 
arate a z = +1 hexagons from a z = — 1 ones. The 
actual path of these loops will depend on the Dq 
but for a given reference dimerization only one such 
path exists. The reason for this is easily under- 
stood by looking at a single hexagon: whatever Dq 
maybe there is always a single loop which surround 
this hexagon and this hexagon only. 

We can check the above property by a direct counting. 
On the one hand we have 2 W / 3 pseudospins configurations 
and 2 Ar / 3 ~ 1 pairs of non equivalent configurations. On 
the other hand there are 2 Ar / 3+1 (see Ref.[l|) dimerizations 
on a kagome lattice with N sites and periodic boundary 
conditions. The agreement is found by remarking that 
the number of dimerizations has to be divided by four to 
get the size of a single topological sector 



2. Pseudo-spin flip operator a x 

One interest of ZE pseudospin representation is that 
the er^ operator which flips the pseudospin at position h 
can be expressed in a simple way in terms of local dimer 
operators. It seems that ZE did not realize this very use- 
ful property of their representation. The simplest dimer 
moves involve loops around hexagons. These 32 loops are 
represented in Tabled The corresponding operators 



L a — \L a ) {L a \ + (L a \ 



(Al) 



shift the dimers along the loop L a if it is possible and 
annihilates the state otherwise. We will now prove that 
a x is the sum of all the 32 kinetic operators of hexagon 
h: 



:S2 



(A2) 



a=l 



The fact that this sum of dimer operators L a real- 
izes the spin algebra is not obvious. In particular the 
fact that these operators commute at two neighboring 
hexagons h and h! must be verified since in general 



L a (h),L a ,(h>) 



*0. 



Consider an arbitrary dimerization \D) in the vicinity 
of an hexagon h. The crucial point is that all the kinetic 
operators L a (h) but one annihilate \D). This is a prop- 
erty of the kagome lattice that we already used before: 
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for any given dimcrization one and only one loop can 
surround hexagon h. So \D') = a x (h) \D) is a dimer con- 
figuration which differs from |D) by a single loop around 
h. Using the o~ z base to represent dimer coverings we 
know that such a state is unique and is the state ob- 
tained from \D) by flipping the pseudospin at in h. Thus 
we have shown that a x (h)a z (h) = —<j z (h)a x (h). 

One can use a very similar reasoning to show that 
a x (h) and a x (h') commute when h^h! but this result 
is most easily obtained by the arrow representation. 



3. Counting dimer coverings with pseudospins 

The Ising basis of ZE's pseudospins provides a way of 
counting the number Af of dimer coverings on any lattice 
K which is the medial lattice of a trivalent one. The 
result is 

A/Y • — n2 N P e ~ 1 fATl 

J v dim. coverings — v-^ / 

where N ps is the number of pseudospins and a is the 
number of topological sectors, it is related to the genus 
by a = 2 2 - 9 in the two-dimensional cases. The factor 
2 -1 comes from the fact that a dimer configuration cor- 
responds to two pseudospin states. Using the Euler rela- 
tion we can check that Eo. I A3I indeed coincides with the 
result obtained with the arrow representation (Eq. ^) : 

\f,. - o2 S9 A73+2-2 3 -l _ 9 iV/3+l / \a\ 

J v dim. coverings — ^ z, — V / 



APPENDIX B: fj, KINETIC OPERATORS AND ZE 
PSEUDOSPINS 

In order to write /x with o~ z and a x only we need to 
express the sign e(h) of Eas.ll0landll2lin terms of the a z 
operators on the neighboring hexagons. We will do this 
with the help of the arrow representation. The argument 
generalizes easily to the fl operators. 

First draw the arrow representation of the reference 
dimerization \Dq) in the vicinity of hexagon h, as in 
Fig. Each neighboring hexagon of h has two ar- 

rows which belongs to the star of h. These arrows can 
cither: A\) point toward the exterior of h; A%) point 
toward the interior of h; or B) point in two different di- 
rections. In Fig. 1121 for example, we have no hexagon in 
case A\, 2 in case A^ and 4 in case B. First look at the 
change in n ou t(h) when a single neighboring pseudospin 
is flipped with respect to the reference configuration (it 
has a z = — 1). If the corresponding hexagon is of type 
A\ (rcsp. A2), n ou t(h) is decreased (resp. increased) by 2 
units. If that hexagon is of type B, n on t(h) is unchanged. 
On a state where a single pseudospin is down we have 
therefore showed that: e(h) = e ro f (h) Yl i<£A c{ ^ cr z (i) 
where A ie f(h) is the set of the neighbors of h which are 
of type A\ or Ai in the reference dimcrization. e re f(h) is 
the value of e(h) in the reference state. 



Now look at the value of e(h) when the neighbor- 
ing pseudospins are in an arbitrary state a z (hi) = ±1. 
If two pseudospins are flipped on hexagons which are 
themselves first neighbors, one arrow is flipped twice 
and therefore remains unchanged. This adds to e(h) a 

— 1 factor, which multiplies the single pseudospin factor 
discussed above. This can be seen on the example of 
Fig. If the pseudospin of the hexagon containing the 
sites (2', 3, 3') (type A 2 ) is flipped the sign e(h) changes 
but if hexagon (l',2, 2') (type B) is flipped e(h) is not 
affected. However, if both hexagons are flipped e{h) is 
unchanged too. In addition to the a z (hi) factors com- 
ing from hexagons belonging to A le ((h) we must add a 

— 1 contribution for each pair of consecutive pseudospins 
which are simultaneously in a a z {hi) — —1 state. Such a 
factor is given by: 

6 

JJ/[<r*(ft i ),tr*(ft i+ i)] (Bl) 

i=i 

where / is equal to — 1 when the two pseudospins are both 
down a z {hi) = a z (hj + i) = —1 and = 1 otherwise. The 
product runs over the six neighbors hi of h numerated in 
a cyclic way. / can be explicitly written as a polynomial: 

f[af,al] = 1(1 + of +02-0^2) (B2) 

We eventually have an expression for e{h) as a polynomial 
of ZE pseudospins operators located on neighbors of h: 

6 

e(h) = e iei (h) H a z (i)l[f[a z (hi),a z (h i+1 )] 

i£Arei(h) i=l 

(B3) 

The left-hand side is both a local and reference- 
independent operator. On the right-hand side, the in- 
formation on the arrows in the reference configuration is 
present at several places: in e re f(/i), in the set of sites 
A IB f(h) as well as in the a z operators. It is possible to 
check directly on this expression that the e(h)a x (h) sat- 
isfy the [i algebra. 

APPENDIX C: REMOVING THE EXTENSIVE 
DEGENERACY 

To perform the explicit diagonalization, it is useful to 
remove the extensive degeneracy. Let us denote by C a , 
a = 1, . . . N c , the commuting and independent operators 
introduced in Scc. lVE"5l and c a = ±1 their eigenvalues. 
We can decompose the Hilbcrt space in eigenspaces {c Q }, 
using the projectors 

P£ = i(l±C a ). 

It is then sufficient to work within the reduced space 
where c a are all equal to 1, for example. The reduced 
space is then generated as follows. Consider the state 
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|{s(/i) = l};{Ik}) where all s(h) — 1, and the integrals 
of motion have some values {Ik}- This state exists, is 
unique and it has the property of transforming the action 
n(h) into that of p,(h) for all h. We project this state to 
the eigenspace {c a } = 1, 

\n;{I k })=l[PZ\{s(h) = l};{I k }) 

(V 

and then generate the whole subspace {c a = 1} by the 
action of the monomials in [i 

N ps 

Y[rth) n *\Sl;{I k }) , 

h=l 

with rifi = 0,1. Roughly half of the spins can be elimi- 
nated recursively by the following procedure. For spins 
belonging to the sublattice A, we bring n(h) at right us- 
ing the commutation relations, we transform it in p,(h), 



which is 1 by construction. Some of the spins can be 
eliminated using the integrals of motion And finally, 
one of the spins entering to a bow-tie operator can be 
eliminated by replacing it with the product of the other 
three spins in the bow-tie and the operator T(h). The 
latter operator is brought to the right, transformed in 
T(h) whose value is 1. Similarly, one can construct the 
eigenspace for any values of {c a }. The passage from a 
sector of fixed {c a } to another can be realized by op- 
erators fl(h) belonging to one of the sublattices B, C, D. 
Since these operators commute with the Hamiltonian, the 
spectrum of H.^ will be the same in all the sectors, which 
proves that the global degeneracy of the spectrum is 2 Nc . 

The action of the Hamiltonian on the reduced basis 
is computed using the same procedure as above. The 
fact that we work in a space with dimension divided by 
2 n ps/2 allows us to perform numerical diagonalization for 
relatively large systems, up to N ps = 48. 
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differ from those of Mambrini and Mila who worked in the 
full dimer space. 

55 We repeat the fit when one and when any two data points 
are removed. The minimum and maximum values of the 
coefficient are used to estimate the uncertainty. 

56 This move should not, however, contain topologically non- 
trivial loops. 

57 A dimer on bond ij oriented from i to j then corresponds 

to (ITU,) - \hU))/V2. 

58 Notice the difference with bi-partite lattices where it is 
possible to orient the bonds so that the sign (a\b) is always 



be positive—. 

59 The length L max of longest loop around a star is always 
even in our family of lattices since L max = 2L m i n . For this 
reason Ea. 1121 is valid for any such lattices. 

60 It is not surprising to find such a three-sublattice order in 
this variational approximation since the translation invari- 
ance is explicitly broken by Ea. 1341 

61 They are not all independent on a finite system with peri- 
odic boundary conditions, in the sense that some products 
on closed straight lines reduce to products of the integrals 
of motion (monomials in the fi operators - the 1^ defined 
below). The number of dependent operators goes as the 
number of lines M ~ y N ps , so in the thermodynamic 
limit they represent negligible fraction of the total number 
of commuting operators. The degeneracy of the spectrum 
will be 2 N ", of the order of 2 N ^' 2 

62 For Eq. l5Ul to hold, the two bonds must be at a large enough 
distance so that the jl and T involved in the Oi of one 
triangle do not touch the second triangle. This condition is 
automatically satisfied beyond 6 kagome-lattice spacings. 

63 We do not consider the special case when one of the lines 
has length 2. 

64 Two operators and F r anticommute if they correspond 
to two independent cycles, so we can choose F^ and r r to 
generate the four topological sectors on the torus. 

65 It turns out that, as a finite size effect, the N ps = 12 and 
16 systems also exactly have < n >= — i in their ground- 
state. 

06 3 or 6 depending whether the axis symmetry is also spon- 
taneously broken. 

67 If one first neglects the hard core repulsion of dimers, one 
can build a state where all the dimers are in the same 
zero-momentum wave-function <j>{r) = 1. The many-dimer 
wave function is of course ip(ri, ■ ■ ■ , r n ) = <j>(ri) ■ ■ ■ 4>{r n ) = 
1 (already symmetrized) whatever the dimer positions 
ri , • • • , r n are. Projecting this state (Gutzwiller like) onto 
the space of no double-occupancy precisely gives the equal- 
amplitude superposition of all fully-packed dimer configu- 
rations. 

08 It is a property of the kagome lattice that the four topo- 
logical sectors have exactly the same dimension. 



